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Abstract 

W-based granular metals have been prepared by electron beam induced deposition 
from the tungsten-hexacarbonyl W(C0)6 precursor. In situ electrical conductivity mea- 
surements have been performed to monitor the growth process and to investigate the 
behavior of the deposit under electron beam post irradiation and by exposure to air. 
During the first part of the growth process, the electrical conductivity grows non-linearly, 
independent of the electron beam parameters. This behavior is interpreted as the result 
of the increase of the W-particles diameter. Once the growth process is terminated, the 
electrical conductivity decreases with the logarithm of time, a ~ ln(t). Temperature- 
dependent conductivity measurements of the deposits reveal that the electrical transport 
takes place by means of electron tunneling either between W-metal grains or between 
grains and trap sites in the matrix. After venting the electron microscope the electrical 



1 



conductivity of the deposits shows a degradation behavior, which depends on the compo- 
sition. Electron post-irradiation increases the electrical conductivity of the deposits. 

1. Introduction 

Electron beam induced deposition (EBID) is a high resolution one-step technique used 
to deposit and to pattern two- and three-dimensional micro- and nano-structures [1]. The 
importance of EBID is rapidly increasing in applied science and fundamental research [2] . 
On the one hand, the possibility of direct writing makes EBID a promising alternative to 
nano-lithography and a useful tool for mask repair. On the other hand, the capability to 
produce deposits from many different precursors with tunable electrical properties makes 
this technique attractive for the development of new materials. EBID is based on the 
interaction of an electron beam with a substrate which is covered by adsorbed precursor 
molecules and which contain the metal or semiconductor to be deposited. The electrons 
dissociate the precursor molecules into a volatile component, which leaves the surface and 
into a non-volatile one, which forms the deposit. The deposits consist of a disordered 
array of crystalline metallic nanoparticles with diameters between about 1 nm to 5 nm 
embedded in an insulating matrix. The metal volume fraction, i.e. the average particle 
size and the interparticle distance, can be varied by tuning the electron beam parameters 
(beam current, acceleration voltage, dwell time). 

In situ electrical conductivity measurements of granular materials prepared by EBID 
or IBID (ion beam induced deposition) are valuable in order to study the growth pro- 
cess and the response of the material during post-irradiation and exposure to air. By 
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means of these measurements information about the electrical transport properties, the 
microstructure, the chemical and physical stability of the deposit are deduced. In lit- 
erature, in situ electrical measurements of IBID and EBID deposits are rare [2]. The 
electrical behavior of EBID deposits from acrylic acid has been monitored by means of 
two-probe measurements [2]. Recently an ageing process has been monitored in platinum- 
based nanostructures [3]. The authors report a continuous decrease of the conductivity 
over a time range of many days. Studies performed by using W(CO) 6 as a precursor are 
known for deposits produced by IBID [4] and EBID [5,6]. The investigation of Hoyle 
et al. [5, 6] represents an important reference for the present study. In their work the au- 
thors prepared structures for beam energies between 2 keV and 20 keV. They investigated 
in situ the electrical conductivity of the deposits and, by means of transmission electron 
microscopy (TEM), their microstructure. It is the purpose of the present paper to further 
investigate the properties of W-based deposits from the W(CO)6 precursor by performing 
in situ transient electrical conductivity measurements during deposition, post-irradiation 
and exposure to air. 

2. Experimental 

To prepare our samples we used a dual beam SEM/FIB microscope (FEI, Nova 
Nanolab 600) with Schottky electron emitter and an ultimate resolution of 1 nm. In 
this system the electron beam power can be continuously tuned by means of the contin- 
uous variation of the beam energy and pre-defined discrete values of the beam current. 
The microscope is equipped with a gas injection module which introduces the W(CO)6 gas 
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precursor via a 0.5 mm diameter capillary in close proximity to the focus of the electron or 
ion beam on the substrate surface. EBID structures were grown on a Si (p-doped)/Si02 
(300 nm) substrate. The substrates were pre-patterned with 120 nm thick Au/Cr con- 
tacts defined by UV-photolithography. In Fig. [T] we show a scanning electron microscope 
(SEM) image of three deposits for two-probe electrical measurements. This technique was 
chosen after having verified that the influence of the contact resistance between electrodes 
and deposits is below about 3%. Our chips are prepared to allow measurements for up 
to 12 deposits. After deposition in situ energy dispersive x-rays analysis (EDX) at 5 keV 
electron beam energy was performed in order to determine the material composition of 
the deposit. The low beam energy was chosen to avoid excitation of x-ray fluorescence in 
the substrate material. This was verified by Monte Carlo simulations of the electron tra- 
jectories for the given thicknesses and compositions of the deposit [7]. For in situ transient 
electrical conductivity measurements a Keithley Sourcemeter coupled with a Multiplexer 
was used to perform current measurements at fixed bias voltage. The conductivity was 
deduced from the known dimensions of the deposits. The lengths I, the widths w and 
the hight h of the samples were determined by direct analysis of the SEM images. The 
maximum geometry-dependent error for the conductivity data amounts to about 35%. 
Finally, temperature-dependent measurements of the electrical conductivity where per- 
formed in a variable-temperature insert mounted in a 4 He cryostat in the temperature 
range 1.8-265 K. 

3. Measurements 
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3.1 Deposition parameters and EDX characterization 

In Fig. [2] we report the results of the composition analysis performed by means of 
EDX for the samples used in this work. The deposits have been obtained with a variable 
voltage and beam current in the range 4 keV< E <20 keV and 0.25 nA< I <6.6 nA, 
respectively. The corresponding beam power varied between 5 nA-keV< p <26 nA-keV. 
Within this range the deposits are granular with a W-content which increases linearly 
between 8.1 at% and 38.7 at%. In the left inset of Fig. [2] we plot the conductivity vs. the 
dose per scan. For high doses per scan the conductivity lies in the range between ca. 
9300 to 18000 f2 _1 m _1 . These values are in agreement with the ones reported by Hoyle 
et al. [5,6]. For smaller beam power we measure a strong decrease of the conductivity 
which reaches 17 f2 _1 m -1 at 65 C/m 2 . This value is one order of magnitude smaller than 
the one reported in Ref. [5,6] for comparable dose per scan. On the right side of Fig. [5] 
we plot the ratio between oxygen (or carbon) and tungsten ([0]/[W] or [C]/[W]) content 
vs. the tungsten-content for each deposit. The ratio ([0] + [C])/[W] is also plotted. For 
W content equal to 8.1 at% we obtain ([0] + [C])/[C]=11.4. This value is close to the 
one deducible from the stoichiometry of W(CO)6 if one were to assume that the CO 
molecules do not dissociate. By increasing the beam power the W content increases, 
whereas the carbon and oxygen content decrease. Additional information can be deduced 
from the ratio [C]/[0] (see inset). For low doses we find [C]/[0]>1, which shows that the 
composition of the matrix is dominated by carbon. The [C]/[0] ratio shows a peak for a 
W content of ca. 15 at% and then it decreases towards 1. By means of TEM measurements 
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Hoyle et al. [5] found an amorphous material for deposits prepared with a dose per scan 
smaller than 500 C/m 2 . For higher doses they found W containing nanocrystals consistent 
with the high temperature /3-phase of tungsten carbide (/?-WCi_ x ). The nanocrystal size 
was estimated to be less than about 3 nm. From the EDX analysis of our deposits we 
point out that the [W] content increases with increasing beam power or with increasing 
dose per scan. 

3.2 Measurements during growth 

In Fig. [3] we show transient electrical conductivity measurements during the growth of 
deposits with a W content of 8.1 at% (sample #1), 14.7 at% (sample #4) and 38.7 at% 
(sample #9). For details concerning the composition, the geometry and the beam param- 
eters used during deposition see Tab. [TJ 

At the time t = 10 sec the gas precursor enters the vacuum chamber. The e-beam is 
rastered over a rectangular area defining the later deposit. In order to prepare samples 
#1 and #4 we used a dwell time of 100 /is per pixel and a pitch of 20 nm. For sample #9 
we used a dwell time of 5 /is per pixel and a pitch of 20 nm. The conductivity changes 
at a rate a'— cr/t which increases monotonically with time during the deposition of all 
samples. Sample #9 shows a rapid increase of cr' in the first seconds of deposition, which 
we attribute to the formation of the first layer of the deposit. In the first 5 seconds cr' 
follows a power law of the form cr' = t a , with a = 3. After that, cr' drastically decreases 
tending to become constant, as can be expected, if the thickness of the deposit grows 
linearly with time, cr' for the deposits with 8.1 at% and 14.7 at% metal content follows a 
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Table 1: Composition and deposition parameters of the samples. Sample geometry 
(l//xm x w//im x h/nm), beam parameters (dwell time//is, pitch/nm): Sample #1 
(8.75 x 1.03 x 240), (100, 20); #2 (9.19 x 1.03 x 315), (100, 20); #3 (8.27 x 1.06 x 226), 
(100, 20); #4 (4.76 x 1.09 x 228), (100, 20); #5 (6.15 x 1.03 x 148), (100, 20); #6 
(18.1 x 1 x 352), (100, 20); #7 (4.57 x 1.03 x 241), (100, 20); #8 (14.93 x 1.63 x 192), 
(100, 20); #9 (7.75 x 1.51 x 682), (5, 20). 



Sample 


W [at%\ 


C [at%] 


[at%\ 


Energy /keV 


Current / nA 


#1 


8.1 


66.4 


25.5 


20 


0.25 


#2 


8.9 


69.8 


21.3 


24 


0.28 


#3 


9.0 


64.9 


26.1 


20 


0.25 


#4 


14.7 


67.5 


17.9 


20 


0.5 


#5 


15.6 


68.9 


15.6 


20 


0.51 


#6 


16.7 


66.7 


16.6 


11 


2.3 


#7 


31.5 


21.1 


47.4 


17 


1.52 


#8 


31.8 


44.4 


23.8 


5 


3.7 


#9 


38.7 


34.7 


26.6 


4 


6.6 



power law of the form <r' ~ t a , with a = 1.86 and a = 1.43, respectively. For sample #1 
cr' follows this power law during the whole deposition process, for the sample #4 it shows 
a decrease after about t = 10 3 s. The fit to sample #4 was made between 300 and 1000 
seconds. In this range the growth rate is lower than in the first 300 seconds of deposition, 
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which we excluded from the fit because of the high dispersion of the experimental points. 
Therefore it is not surprising that the result of the fit gives for a a value larger for sample 
#1 than for sample $4, which is counterintuitive since we expect higher growth rate for 
deposits with higher metal content. 

In conclusion, these results clearly indicate that a linear increase of the thickness is not 
sufficient to explain the transient electrical conductivity measurements during the growth 
of the deposits. Since the increase of the conductivity is faster than linear, an additional 
mechanism has to be considered in order to understand this observation. Among the 
possible mechanisms which may contribute to the enhancement of er' we consider the 
charging effect and the increase of the metal grain size due to the electron beam irradiation. 
In the last chapter of the paper we discuss the relevance of these mechanisms for the 
samples prepared in this work. 

The growth process terminates with the simultaneous shut-off of the electron beam 
and the precursor gas supply. Correspondingly, the conductivity starts to decrease (see 
Fig. [3]). The relaxation of the conductivity can be attributed to the migration of excess 
electrons injected by the electron beam towards the electrodes. In Fig. Hh we depict the 
relaxation for samples #1, #4 and #9. The conductivity relaxation follows a logarithmic 
time dependence, i.e., a = b ■ ln(t). With suitable choice of the parameter b, this 
formula can be used to fit the relaxation of all the deposits. The velocity of the relaxation 
depends on the deposits' composition. In particular, we notice that deposits with lower 
W content relax more rapidly than deposits with higher W content. This information is 
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summed up in Fig. [lb for all the samples prepared in this work. 

3.3 Measurements during venting the microscope 

Transient electrical conductivity measurements during exposure to air where per- 
formed after deposition by venting the electron microscope. The result of these experi- 
ments are summarized in Fig. [5] where we plot the normalized conductivity vs. time for 
various deposits. The venting procedure starts at t=60 s, when iV 2 gas enters the cham- 
ber. The conductivity of the deposits starts to decrease with a rate which depends on 
the metal concentration. In particular, the lower is the metal content the faster is the 
decrease. The reduction of the conductivity becomes faster at £^400 s, when the door of 
the chamber slightly opens and a small flux of air enters the microscope. At fc«900 s the 
door of the microscope is deliberately opened to fully expose the deposits to air. At this 
point the reduction rate of the conductivity strongly increases. According to Fig. 0, we 
divide the deposits in two groups exhibiting strong and weak degradation of the conduc- 
tivity, respectively. To the first group belong the samples #1, #2, #5 with metal content 
equal to 8.1 at%, 8.9 at% and 15.6 at%. After one hour from opening of the vacuum 
chamber the conductivity of these deposits has dropped to 7 % to 34 % of its initial value. 
To the second group belong samples #6, #7, #9 with 16.7 at%, 31.5 at% and 38.7 at%. 
In this case after one hour of exposure the value of the conductivity lies between 90 % 
and 99 % of the initial value. It is interesting to note that the degradation rate increases 
monotonically with decreasing metal content and, thus, with the beam power used to 
prepare the respective samples. Most likely the degradation is due to a reduction of the 
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tunneling probability within the deposit. In the discussion section we speculate that the 
degradation rate can be linked to the density of the deposit which depends on the beam 
power [8]. 

3.4 Post- irradiation conductivity behavior 

In order to study the effect of the beam irradiation on EBID deposits we consider 
two samples with W content of 9 at% (sample #3) and 14.7 at% (sample #4), respec- 
tively. After preparation the samples were irradiated with the same beam power used 
during deposition, i.e. 5 nA-keV and 10 nA-keV, respectively. In Fig. Eh we plot the tran- 
sient conductivity vs. the irradiation time. The conductivity of both samples abruptly 
increases during the first few seconds of irradiation. In the following minutes the conduc- 
tivity increases monotonically at a much lower rate, which depends on the composition of 
the deposit. The conductivity of sample #3 increases to approximatively eight times its 
initial value in 2000 seconds of irradiation. In the same time interval the conductivity of 
sample #4 increase only by 10 % of its initial value. The increase of conductivity may be 
attributed to the increase of charge carriers released during the beam-induced breakage 
of carbon-carbon bonds in the amorphous matrix, as was speculated in Ref. [9, 10]. Sam- 
ple #3 shows the largest increase of conductivity because it is the sample with largest 
matrix volume fraction. In Fig. [6b we report the relaxation of the conductivity over a 
time scale of almost 3 days. The most remarkable fact is the very low decrease of the con- 
ductivity in comparison with its increase during post-irradiation. The second remarkable 
fact is that the conductivity is not logarithmically dependent on time. Therefore no relax- 



10 



ation of the conductivity takes place; rather, we suggest that the conductivity deteriorates 
because of the exposure to the rest gases present in the SEM chamber. The phenomena 
is similar to the one occurring by exposure to air. However, its time scale is much longer 
because the residual pressure (~ 4-10 _6 mbar) inside the microscope is order of magnitude 
lower than in the environment. Note that in the frame of the "bond-breaking" -model we 
expect that carbon atoms bond again leading to a conductivity decrease much faster than 
the one observed. Therefore we attribute the increasing conductivity during irradiation 
to a permanent chemical transformation of the deposit. 

3.5 Temperature dependence of the electrical conductivity 
In this section we present the temperature-dependent conductivity measurements per- 
formed for the samples with W content of 8.1 at% (sample #1), 14.7 at% (sample #4) 
and 31.8 at% (sample #8). For the low conductivity samples #1 and #4 we choose an 
applied field between 1143 V/cm and 210 V/cm in order to be able to measure cx(T) 
down to the lowest temperature accessible in our setup. In order to avoid irreversible 
switching in the I(V) characteristic during the measurements we used a current limit. For 
sample #8 we choose an applied field of 6.7 V/cm in order to measure cx(T) in the linear 
regime of the I(V) characteristic. In Fig. [7^ we plot the conductivity vs. temperature 
dependence for samples #1 and #4. Both samples show the characteristic of systems 
governed by variable-range-hopping, i.e. a = cr exp (T /T) a . In particular, we observed 
a = 0.28 and a = 0.37 for samples #1 and #4, respectively. The behavior in the system 
with lowest metal content is close to that one of amorphous carbon films [11, 12]. In that 
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case the tunneling takes place between localized sites in the limit of negligible Coulomb 
blockade effects (a = 0.25) [11]. The value of a obtained for sample #4 may suggest 
both, a contribution to the tunneling between the trap sites in the matrix and between 
the metal grains; the latter mechanism is described by the exponent a = 0.5 [13, 14]. In 
Fig. [7b we plot the temperature dependence of the electrical conductivity for sample #8. 
We find that the conductivity follows a power law over the complete temperature range 
of our measurement (1.8 K< T <265 K) with a = a Q + a T@, with f3 = 0.55. Such a 
behavior is similar to the one recently reported in Ref. [15] and interpreted by means of 
a tunneling percolation model in the limit of large inter-grain tunneling [15, 16]. 

4. Discussion 

A linear increase of the time-dependent conductivity can be explained by supposing 
a linear time dependent growth of the thickness. However, the measurements performed 
during growth show an increase of the conductivity vs. time faster than linear. Therefore, 
an additional mechanism which contributes to the enhancement of the conductivity has to 
be considered. Possible additional mechanisms for the increase of the conductivity growth- 
rate include an enhancement of the electron carrier density and of the metal-particles' size 
during deposition, both due to the electron beam irradiation. In the following we relate 
the time dependence of the conductivity to the increase of the metal particles' size by 
means of a "shell-model". We consider a sphere of radius r and volume v s representing 
a metal-particle, see Fig. [HI Coordination sphere of radius I and volume v i surrounds the 
first one. The radius difference s=l-r represents the separation between two neighboring 
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particles. The ratio v = v s jv\ gives the metal- volume fraction of the system. It follows 
that the diameter of the metal particle is d = 2 r = 2 I i>5. We suppose that at time 
t=to the deposit is constituted by particles with the same radius r. The particles of this 
initial deposit increase their volume because of electron beam stimulated diffusion of W- 
atoms. In general, the growth rate of the particle size is described by the Lifshitz-Slyozov 
coarsening law r ~ £3 [17]. However, also lower growth rates with r ~ i« with n>3, 
have been reported [18,19], which indicate diffusion limited mass transport. Known the 
diameter d, the separation s and the growth rate, the time-dependence conductivity in 
the VRH regime can be deduced from equation (pQ) 

a ~ e * k B T (1) 

where £ is the localization length, ks the Boltzmann constant, Tthe temperature and 
W the effective activation energy for charge transport [20,21]. The time dependence of 
the conductivity is implicit in the localization length and in the activation energy, both 
dependent on the particle diameter. In particular, £ = ^ mx {l + d/s), with ^ mx being the 
localization length of the matrix [20] . The localization length is obtained from the relation 
^ mx = 2.8e 2 /4ne€ok B T (l + d/s) [22], where e is the dielectric constant of the deposit and 
T a constant deducible from the temperature-dependent conductivity measurements. 
From the experiments one also extracts the activation energy W=a kT(T Q /T) a [20]. From 
FigEK it results T = 1.28 TO 5 K and T = 6.4 -10 3 K, for samples #1 and #4, respectively. 
The values of the dielectric constant necessary to calculate £ mx are chosen considering that 
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for a carbon-matrix e ~4 [23] and that with increasing metal volume fraction, the dielectric 
constant of a metal- insulator composite increases [24]. By using for e the values 8 and 
12, it follows £ mx = 0.02 nm and £ mx = 0.1 nm, £=0.03 nm and £=0.61 nm for samples 
#1 and #4, respectively. For a comparison note that the value of the localization length 
of amorphous carbon thin films lies between 0.1 nm and 1.2 nm [25]. In Fig. [9] we report 
the time-dependent conductivity measurements for samples #1 and #4. In order to fit 
the experimental data, we choose a linear increase of the metal volume fraction v from 
0.04 to 0.06 and from 0.04 to 0.35 for samples #1 and #4, respectively. The best fit 
was obtained with n=9.8 and n=6.3, 0.61< s/1 <0.65 and 0.28< s/1 <0.66, respectively 
for samples #1 and #4. These choices for / and v lead to 0.55 nm < s < 0.59 nm 
and 0.42 nm < s < 0.99 nm, 0.62 nm < d <0.71 nm and 1.03 nm < d < 2.16 nm, 
respectively for the two samples. In conclusion, the experimental data on the growth 
process can be well described within our " shell-model" by using a coarsening growth rate 
which depends on the deposit's metal content and, thus, on the beam parameters. We 
find that in the deposit with lowest metal content (sample #1, n=9.8) the size of the 
metal particle increases in time with the lowest rate, as one may have expected. The 
coarsening law of Lifschitz-Slyozov is not reached by the range of beam parameters used 
for sample #1 and ^4. We believe that for higher deposition rates, such regime may be 
reached. It has to be noted that we cannot describe the behavior of sample #9 by means 
of the previous modelling because the parameters T and £ cannot be deduced from the 
temperature dependence of the conductivity of Fig. [7b. However, it is quite possible that 
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in the first stages of the deposition process the electron transport may be described by 
the VRH mechanism. In this case, the growth rate of the particle size of sample #9 shall 
be much higher than the one of samples #1 and $4, being close to the Lifschitz-Slyozov 
theory. Indeed, an hint that this may be the case is given by the different power law, 
which is a= t 4 for sample #9 in the first seconds of deposition, while it is a= t 2 86 for 
sample #1, during the whole range of deposition time. It is important to note that the 
"shell-model" assumes the presence of nanocrystal particles in the deposit. According to 
the TEM measurements of Hoyle et al. [5] this assumption is fully justified for deposits 
prepared with high doses. For deposits prepared with a dose per scan smaller than 500 
C/m 2 the authors of Ref. [5] found amorphous material. However, the presence of very 
small crystals not detectable from TEM measurements cannot be excluded. Indeed, the 
sharp fall of the conductivity with the decrease of the doses (see inset in Fig. [2]) may also 
be explained by the Coulomb blockade of small nanocrystals at room temperature. 

In granular materials it is observed that in the low conductivity tunneling regime the 
temperature dependence of the electrical conductivity is of the form a ~ exp (To/T) a , 
with 0.25< a <0.5. For our samples we find 0.36< a <0.5 with metal content < 19%, 
which may indicate a tunneling mechanism either involving localized sites in the matrix 
or directly between the metal grains. However, in order to have tunneling the typical 
hopping distance for the electron, r* = ( Awe e eo i T ) ' 5 , has to be larger than the distance s 
between the localized states [22]. For the samples #1 and $4 we obtain r*=0.46 nm and 
r*=1.69 nm, respectively. By comparing these numbers with the ones for s given above 
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by the "shell- model" we find that in our deposits direct tunneling between grains appears 
to be possible. On the other hand, tunneling between localized sites in the matrix is not 
favorable. However, this does not exclude conduction inside the matrix which may also 
take place between localized matrix sites bridged by metal grains. As a consequence in 
the following we speculate that the transport mechanism in our deposits takes place by 
means of two channels: directly between grains and indirectly between localized sites in 
the matrix through a grain. In order to describe the transport mechanism we use the 
following formula: 

a = a ■ ooi ■ e ( t > + b ■ a m • e 1 T ' + a 03 ■ e ( t > (2) 

The first two terms of equation (j2J) take into account the tunneling of electrons between 
trap sites in the matrix and between grains, respectively. The factor a gives the volume 
concentration of the matrix, while b refers to the concentration of the metal. The third 
term of the equation takes into account the hopping of electrons between trap sites and 
grains. The temperature T03 and factor 5 are free parameters. By fitting the experimental 
curves with equation (J2]), we obtain the volume concentration of the matrix and of the 
metal. We made the fit to the two curves of Fig. [7] ('"fit 1"' and '"fit 2"'), for which we 
have W content of 8.1at% and 14.7at%, associated to a = 0.39 and a = 0.36, respectively. 
For the fit we set a 01 = 10 Q ~ l cm" 1 , T 01 = 10 7 K [25], a 02 = 27 Q ~ x cm' 1 and T 02 = 
2070 K [15]. For the third term of eq. © we choose a 3 = 15 Q ~ l cm' 1 , T 03 = 1.3-10 4 K 
(fit 1) and T 03 = 7-10 3 K (fit 2). The best fit was obtained for a = 0.06, b = 0.94 and 
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S = 0.38 (fit 1) and a = 0.35, b = 0.65 and 5 = 0.38 (fit 2). We can compare the 
fractional volume obtained by the fit with the one deduced by EDX measurements. The 
fractional volume occupied by the metal particles depends on the density of the metal 
and of the matrix, which is not well known. The lower and upper limit corresponding to 
a density of 1.8 g/cm 3 (graphite-like amorphous) and of 3.5 g/cm 3 (diamond-like) lead 
for the W-content of 8.1at% to a fractional volume between 12vol% and 20vol% [15]. 
From our fit we obtain the volume of 6%, which would indicate a lower density than 
a graphite-like matrix. This difference may be due to the presence of oxygen atoms 
inside the carbonaceous matrix. By means of in situ mass measurements the density of 
the matrix of a W-based deposits from W(CO) 6 has been related to the electron beam 
parameters used during deposition [8]. In general, the density of the matrix increases 
with the power of the electron beam. The authors measured a matrix density between 
0.29 g/cm 3 and 0.88 g/cm 3 corresponding to a W- volume concentration between 11.1 vol% 
and 6.4 vol% [8,26]. These results where obtained for beam voltages between 5 keV and 
15 kev, but with a beam current much higher than the one used in our depositions. For 
a system with W-content of 14.7at%, we calculate a fractional volume between 19vol% 
and 33vol% for graphite- and diamond-like matrix, respectively. The value of 35 vol% 
obtained within our modelling is close to that for a diamond-like matrix. Of course the 
matrix of our deposit is non-crystalline and does also contain oxygen. In conclusion, we 
find clear evidence of an increase of the density of the matrix by increasing the electron 
beam power. This result is in accord with the one obtained for amorphous carbon deposits 
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obtained by electron-beam-induced deposition [8,26]. 

After terminating the growth process the conductivity starts to decrease with the log- 
arithm of the time, see Fig. HI The relaxation can be attributed to the migration of the 
charge carriers out of the sample. In general, in order to describe the relaxation of a 
system toward equilibrium a stretched exponential of the form a ~ exp (— t/r)' 3 , with 
0< (3 <1 can be used [27]. In the case of 13=1 the the Debye relaxation is obtained, which 
is valid for systems with weak dynamic correlations. In the case of hopping conductiv- 
ity, this corresponds to a system with no Coulomb interactions. If the system presents 
strong dynamic interactions, the law of Kohlrausch [27] can be obtained by considering a 
serial relaxation of hierarchies of correlated systems [28] . In case of hopping conductivity, 
this corresponds to a system with strong Coulomb interactions, dynamically described by 
the tunneling of an electron into a localized state, which is possible if a second electron 
tunnels into another localized state (co-tunneling). The logarithmic decay of the variable 
describing a system finds an explanation in the frame of hierarchically constrained dy- 
namics [29]. In particular, in the theoretical study of Ref. [30] it is mentioned that the 
conductivity of systems governed by hopping conduction can exhibit a logarithmic time 
dependence, as we find experimentally. Finally, we observed an increase of the relaxation 
rate at lower W-concentration. We attribute this behavior to the increased number of 
localized sites in the matrix and, therefore, to an higher probability for the electron to 
tunnel into a free localized state. 

The increase of the density of the matrix with the electron beam current (see above) 
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can be used to interpret the conductivity measurements obtained during air exposure. 
The degradation rate a' decreases as one considers deposits prepared at larger beam 
current (power) and, thus, with the density of the matrix. In order to better describe 
the degradation process, we fit the experimental data of Fig. El from the time where the 
door of the microscope is opened (i~900 s), with a sum of two exponential functions: a oc 
ai exp(t/Ti)+a 2 exp(t/r 2 ). In Fig.fTUa we plot a typical decay fitted with such a curve. In 
Fig. [TUb we report the decay constants T\ and r 2 for the deposit measurements depicted in 
Fig. [5] as a function of the beam power used to prepare the samples. Both decay constants 
increase with the beam power with (15 < T\ < 236) s _1 and (510 < r 2 < 2185) s -1 . As one 
observes from the inset of Fig. [TUb the weight of the first term of the exponential sum is 
larger than the weight of the second one (ai>a 2 ). In particular the ratio strongly increases 
for small values of the beam power, which corresponds to deposits mainly constituted by 
the matrix. Therefore we associate the first term of the exponential sum to the matrix. 
The decay constant T\ is associated to the decrease of the tunneling probability within the 
matrix, or between matrix and grains in view of our previous discussion. The increase of 
the decay constant T\ with beam power is interpreted as due to the increase of the matrix 
density. The second term of the exponential sum may describe the decay of the tunneling 
probability between metal grains. The increase of the decay constant r 2 with beam power 
may be interpreted as due to the increase of the grain size. 

Similar considerations are useful in order to analyze the behavior of the deposits af- 
ter irradiation. Also in this case the conductivity can be described by the sum of two 
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exponential functions. The two decays are independent processes linked to the decrease 
of the tunneling probability within the matrix and between the metal grains. The decay 
constants of sample #4 and sample #3 are Ti=1.91T0 4 s, ti=2.33-10 4 s and r 2 =1.66-10 5 s, 
r 2 =1.6-10 5 s, respectively. These values are 2 to 3 orders of magnitude larger than those 
obtained by venting the system. From the literature it is known that the presence of 
residual oxygen and water molecules inside a SEM is sufficient to oxidize W-based gran- 
ular metals prepared from the tungsten hexafluoride (WF 6 ) precursor [31]. Similarly, the 
oxidation of the deposits is most likely the reason for the observed decrease of the con- 
ductivity that we measured inside the SEM. As can be expected, the degradation takes 
place on a much longer time scale in comparison to the exposure to air. 

5. Conclusions 

Probably, the most relevant technological aspect concerning our study is the investi- 
gation of the behavior of W-based granular metals during exposure to air. We have found 
that the degradation of the system strongly depends on the composition. In particular, 
after one hour of exposure to air the remaining conductivity varies between 7% and 99% 
of its initial value for deposits with W-content in the range 8.1at%<W<38.7at%. Tran- 
sient electrical conductivity measurements indicate that this dependence is due to the 
different densities of the deposits. Furthermore, the decay of the conductivity involves 
the matrix and the W-metal grains on two different time scales. In conclusion, in view of 
the possible employment of EBID W-based deposits for applications the composition of 
the composite has to be carefully chosen to minimize ageing effect. In this direction it is 
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remarkable that after one year from deposition two samples with W-content of 34% and 
19%, respectively, show a decrease of the conductivity only of a factor between two and 
four. However, covering the deposit with a protective SiO x , layer like tetraethylorthosil- 
icate (TEOS) precursor, may be necessary to prevent long time scale degradation of the 
electrical properties [3]. 

By means of transient electrical conductivity measurements we have monitored the 
growth process of W-based deposits. In order to explain the non-linear growth of the 
conductivity vs. time we have used a "shell-model" based on the hypothesis that the 
diameter of the W-metal particles grows during deposition. Within the modelling we have 
estimated the average particle size in the final deposit, which ranges between 0.71 nm and 
2.16 nm for deposits with W content of 8.1at% and 14.7at%. In would be interesting to 
integrate this analysis with transmission electron microscopy investigations. 

Furthermore, we have investigated the electrical transport mechanism as a function of 
the deposits' composition by means of conductivity vs. temperature measurements. From 
the analysis of the experimental data we have speculated that the electrical transport takes 
place by means of two parallel channels involving both the direct tunneling between metal- 
particles and the indirect tunneling between localized sites in the matrix bridged by metal 
particles. As a consequence of this investigation information about the matrix density 
can be deduced, which is useful to interpret the data of the conductivity degradation 
under exposure to air. Finally, we have observed that the relaxation of the charge carriers 
after termination of the deposition process follows a logarithmic time dependence. Such 
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a behavior suggest that at the end of the deposition the electron transport in the deposits 
is dominated by strong Coulomb interaction and that the relaxation can be described in 
the framework of hierarchically constrained dynamics commonly used to describe strongly 
interacting glassy materials. 

In conclusion, we have shown that by means of in situ transient electrical conductiv- 
ity measurements very useful information about the electrical transport mechanism, the 
microstructure and the chemical and physical stability of W-based granular materials can 
be obtained. It shall be interesting to extend this measurements to materials grown with 
other precursor gases and to further these studies with TEM investigations. 
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Figure 1: Upper figure: SEM image of three EBID structures for two-probe electrical 
measurements. The deposits are written on 120 nm thick gold/chromium electrodes. 
Lower figure: Setup for the in situ transient electrical conductivity measurements. On 
this chip 3 from 12 possible measuring positions are contacted via this copper wires 
attached to the contact pads with silver paint. 
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Figure 2: Left figure: Dependence of the W- metal content on the electron beam power. 
The conductivity increases with the dose used during deposition, see inset. Right figure: 
EDX analysis of the composition of the samples prepared in this work. The concentration 
of carbon and oxygen decrease with increasing power of the electron beam. In the inset 
we report the ratio between carbon and oxygen. 
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Figure 3: Electrical conductivity measured during the growth of sample #1 (8.1 at% W), 
#4 (14.7 at%) and #9 (38.7 at%), respectively. In each case, the deposition was started at 
t=10 sec. The data points in the inset for sample #9 show that after an initial non-linear 
growth, the conductivity tends to increase linearly with the time. The conductivity of 
samples #1 and #4 is non-linear during the whole deposition time considered, which is 
due to the lower growth rate of these deposits. The electrical conductivity of each sample 
decreases once the deposition process is terminated. 
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Figure 4: a: Relaxation of the conductivity after the end of the deposition process. The 
conductivity decreases with the logarithm of the time, o = b ■ ln(t). b: Relaxation rate 
vs. metal content. The lower the metal content the faster the relaxation rate. 
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Figure 5: Electrical conductivity vs. time as the SEM is vented. During the procedure 
first nitrogen enters the chamber; after ca. 6 minutes the door of the microscope slightly 
opens. 8 minutes later the door is fully open. The conductivity decreases with time at a 
rate which depends on the material composition. See text for details. 
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Figure 6: Post-irradiation experiment. Left figure: Conductivity vs. irradiation time. The 
increase of the conductivity is the highest for the sample with lowest W content. Right 
figure: evolution of the conductivity after termination of the irradiation. The decay of 
the conductivity is described by two parallel decay channels involving the matrix and the 
W grains. The fastest decay is for the deposit with the lowest W content. See text for 
details. 
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Figure 7: Temperature dependence of the electrical conductivity. Left figure: the tempera- 
ture dependence for samples characterized by variable-range-hopping, a oc exp(T /T) a . 
The exponents a=0.39 and a=0.36 are found for samples with W-content of 8.1 at% (sam- 
ple #1) and 14.7 at% (sample #4), respectively. Right figure: temperature dependence 
for sample #8 with W-content of 31.8 at%. The conductivity increases with temperature 
as a = a + aT 13 , with (3 = 0.55, a =5.27 ft^cm" 1 and a=4.08 JTW^K -1 . 
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Figure 8: "Shell-model" used to model the time-dependent conductivity during the growth 
of a deposit. The model assumes that the metal W-particles have initially a radius r. The 
radius growth by Ar in the time At during the deposition, because of e-beam stimulated 
diffusion of W atoms, s is the separation between neighboring particles, I is the sum of s 
and r. 
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Figure 9: Conductivity vs. time during the growth process. The fit to the experimental 
data are made in the frame of the variable-range-hopping theory, see eq.( [T]). The time 
dependence of the conductivity is implicit in the localization length £ and in the effective 
activation energy W, according to the shell model of Fig. [8] 
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Figure 10: Decay of the conductivity during exposure to air. Left figure: the decay, 
as exemplary shown for sample #7, is described by the sum of two exponential decay 
channels according to a — ai exp(t/T 1 )+a 2 exp(t/r 2 ). Each decay channel is associated 
to a tunneling mechanism: either between grains and trap sites in the matrix (decay 
constant Ti) or between metal grains (decay constant r 2 ). Right figure: value of the decay 
constants vs. beam power. The inset shows the ratio between the weight of the two 
tunneling mechanisms: ai/a 2 . For deposits prepared at small beam power the tunneling 
involving the matrix is dominant. 
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